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Abstract 

Exactly solvable variable parametric Burgers type equations in one-dimension are intro- 
duced, and two different approaches for solving the corresponding initial value problems are 
given. The first one is using the relationship between the variable parametric models and their 
standard counterparts. The second approach is a direct linearization of the variable parametric 
Burgers model to a variable parametric parabolic model via a generalized Cole-Hopf transform. 
Eventually, the problem of finding analytic and exact solutions of the variable parametric mod- 
els reduces to that of solving a corresponding second order linear ODE with time dependent 
coefficients. This makes our results applicable to a wide class of exactly solvable Burgers type 
equations related with the classical Sturm-Liouville problems for the orthogonal polynomials. 

1 Introduction 

The nonlinear diffusion equation was known to Forsyth [TJ, and had been discussed by Bateman 
in connection with various viscous flows [2]. J.M. Burgers considered this equation as a model 
of turbulence [3l 0], Hopf and Cole discussed the Burgers equation in context of gas dynamics 
[5], [5], and Lighthill in acoustics [TJ. Recently, Burgers-KPZ turbulence was addressed in [5] and 
[5], where the Burgers equation was considered also as a good approximation to understand the 
formation and distribution of matter at large scales. 

Mathematically, the standard Burgers equation (BE) V t + VV X = vV XXl is one of the best 
known integrable evolution equations, admitting Lax representation, Painleve property, Backlund 
transformation, etc, [lOj . It is a member of the infinite Burgers hierarchy of evolution equations 
which is linearizable in terms of the heat hierarchy, and analytic solutions can be easily found via 
the Cole-Hopf transform. Various type of exact solutions important in theory and applications 
are known for the BE. For example, when the nonlinear effect and the effects of dissipative nature 
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balance each other, BE exhibits shock and multi-shock solitary wave solutions, which can fuse and 
therefore describe non-elastic interactions [TT], [T^]. On the other hand, it was shown that the 
motion of the poles of rational type solutions of BE corresponds to the motion of classical particles 
interacting via simple two-body potentials. Then, generalization of this idea was used to discover 
new integrable many-body problems, see [13] . [14] . 

Recently, there is an increasing interest in nonlinear evolution equations with variable coeffi- 
cients. In particular, for BE the variable coefficients are able to reflect inhomogeneities in media, 
nonuniformities of boundaries, while the external forces, can model for example time evolution 
of the profile of growing interfaces, [TS]. That makes the forced Burgers models with variable 
coefficients more realistic and significant for applications in various fields. In general, however, 
such models are not integrable and rarely admit exact solutions. For some positive results in this 
direction, one can see [15]. [T7] . p] . [T9] . [20] . [2T] . 

The main subject of this work is the one-dimensional variable parametric Burgers equation of 
the form 

U t + ^rlu + UU x = -^--U xx -uj 2 (t)x, -co < x < oo (1) 
H(t) 2fj,(t) ' J 

where T(t) — jj,(t)//i(t) is the damping term, D(t) = l/2/i(t) is the diffusion coefficient, and 
F(x,t) = —uj 2 (t)x is the forcing term which is linear in the space variable x. 

This model was introduced in our previous work [22] , where exact solutions such as shock and 
multi-shock solitary waves, triangular waves, N-waves and rational type solutions were obtained. 
Here, we give an alternative way for finding these solutions and show the equivalence of both 
approaches. From this point of view, our results can be seen as complimentary to the previous 
ones. On the other side, this work is an extension which includes also the corresponding variable 
parametric potential Burgers and the linear parabolic equation. Thus, we provide a complete 
picture showing the relationships between the linear, the nonlinear, variable parametric and the 
standard models related with BE (JXJ) . 

In Sec. 2, first we show that the parabolic problem with variable coefficients and quadratic 
potential, which appears as a linearization of BE ([TJ , can be reduced to that of solving the classical 
heat equation and a corresponding linear ODE. Second, we solve the IVP for the variable parabolic 
model by finding explicitly the evolution operator as a product of formal exponential operators. We 
remark that, the parabolic models discussed in this article are considered mainly as a linearization 
of the Burgers type models. However, solving some problems with specific initial conditions, we 
show that the results of Sec. 2 can be used also independently to discuss some heat or diffusion 
processes. 

In Sec. 3 and Sec. 4 respectively, we consider the variable parametric potential Burgers and 
Burgers type models. We obtain exact solutions of them using two different approaches. The first 
one, is transforming the variable parametric Burgers type model to a standard Burgers equation, 
which in turn can be linearized in the form of the classical heat equation. This approach was applied 
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in previous work |221. where different type of exact solutions were given and analyzed for models 
with constant damping and exponentially decaying diffusion coefficient. The second approach is 
using a generalized Cole-Hopf transform, so that the variable parametric Burgers model is directly 
linearized in the form of the parabolic equation with time variable parameters discussed in Sec. 2. 
The idea is similar to that which we used to construct exactly solvable Schrodinger-Burgers models 
for complex velocity, which linearization takes the form of a Schrodinger equation for harmonic 
oscillator with time variable parameters, |23j . [24] . Finally, we apply our methods to solve some 
specific IVPs, and illustrate the formal procedure. In Sec. 5, we give a scheme which summarizes 
the main results of the article. 



2 Variable parametric parabolic type equation 

In this section, we consider a variable parametric parabolic model 

f ^ = 1 | Kt)u 2 (t) 2 

dt 2fi(t) dx 2 2 ' (2) 

$(x,t)\ t=to = $(x,t Q ), -o° < x < oo, 

related with the IVP for a second-order linear differential equation of the form 

r + ^-r + u 2 (t)r = 0, r(t )=r ^0, r(t ) = 0. (3) 

In Proposition 12.11 we show that one can solve the IVP ^ in terms of solution to the classical 
heat equation and solution r(t) of the IVP ©. 

Proposition 2.1 The IVP for the variable parametric parabolic equation 0) has solution of the 
form 



lr(t Q ) ^ ^ ( /z(t)r(t) , 



where r(t) satisfies the IVP the auxiliary functions 



are 



V (x,t) = r -^lx; r(t)=r 2 (t ) f f r(t ) = 0, (5) 

r(t) J ti{0r 2 {£) 



and tp(r), r) is solution of the IVP for the classical heat equation 



(6) 



ip(ri, 0) = $(7?, t ), -00 < r) < 00. 

Proof: Using the ansatz = exp[(g(t) — p(t)x 2 )/2] ip(e 9 ^x,r(t)), one can show that, if 

the auxiliary functions p, r, g satisfy the nonlinear system of ordinary differential equations 

P+j^ P 2 +Kt)u 2 (t) -0, p(t Q ) = 0, (7) 
r-^=0, r(i o ) = 0, (8) 
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9 + -7TT = 0, g(t ) 

n(t) 



0. 



(9) 



then the IVP J2]| reduces to IVP (|6|). Note that, Eq.© is a nonlinear Riccati equation which can 
be linearized in the form of Eq.©, using p(t) — fi(t)f(t)/r(t). Then, to find g(t) one needs to solve 
g(t) +r(t)/r(t) = 0, g(t ) = 0. For r(t) ^ 0, it has real-valued solution g(t) = - In \r(t)/r(t Q )\, 
but in this work we shall allow g(t) to be complex- valued, that is g{t) = —\n(r(t)/r(to)), which 
leads to more natural results for the Burgers type equations. Then, the system is easily solved, 
and we obtain the auxiliary functions in terms of solution r(t) to the IVP ([3]), that is 



r(t) 



, r(t o )=0; 



Writing these functions back in the ansatz, gives solution (j4]). 

□ 

Since the IVP (JSJ) for the classical heat equation has solution 



r(t) 



1 



exp 



2t 



(10) 



then according to Proposition [2Tl the solution of the IVP ((2]) is found in the form 

(f](x,t) - yf 



l r(to) 
r(t) 



x exp 



2r(t) ' 



1 



exp 



2r(t) 



$(y,t )dy, (11) 



where r(t) is solution of IVP ([3]), and rj(x,t), r(t) are as defined in ([5]). In case &(x, to) — S(x), 
one obtains the fundamental solution 



$(aj,t) 



r(t ) 



x exp 



2nT(t)r(t) 

Solution (|lip can be written also in the form 



_ ( r 2 (t Q ) f i(t)r(t) 
\2 T (t)r 2 (t) 2r(t) 



K t (x,y)&(y,t )dy, 



where the integral kernel is 
K t (x,y) = 



K*o) 



2-KT{t)r{t) 



x exp 



' 2r(t) X . 



x exp 



{r)(x,t) - yf 



2r(t) 



(12) 



(13) 



(14) 



Using Proposition ^. 11 we solve two linear diffusion problems. 

Heat IVP-A. First we consider the well-known constant coefficient heat equation with quadratic 
potential 

9$ 1 <9 2 $ Cj 2 

(15) 



dt 2 dx 2 
$(i,t)| t=0 = $(aj,0) 



-00 < X < 00, 



where fi(t) — 1, co 2 (t) — —uj 2 , Cj > 0. To find the solution and integral kernel of this problem, 
there are many different methods in the literature, one can see for example Ref. |25j . p. 286, and 
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Rcf.[26 , p. 153. As an alternative, we use the approach in Proposition 12.11 The related IVP is, 
r — uj 2 r = 0, r(0) = tq ^ 0, f(0) = 0, and the auxiliary functions are 

r(t) = rocosh(cl)i), T}{x,t) = sech(ujt)x, r(t) = tanh((Dt)/a). 

Then, according to solution of the IVP ([T5j) is found in terms of solution <p of the classical 
heat equation as 



*(a>,t) 



1 



x exp 



— u> tanh((Di)x 2 



cosh(u)i) 

Using (fTTj) . the general solution takes the form 

1 



x tp sech(Qt)x, 



tanh(d)t) 



ui 



27r sinh(a)i) 



x exp 



exp 



[sec/i(wi)a; — y] 



2 1 



2 tanh(wi) /w 



and in case ^(x, 0) = 5(x), one obtains the fundamental solution 

„2 



27r sinh(d»i) 

Solution (IT6l). can be written also in the form 



x exp 



u> x 



2 tanh(wt) 



$(2/,0)dy, (16) 



(17) 



K t (x,y)$(y,0)dy, 



where 



K t (x,y) 



2n smh(dit) 



x exp 



Q_ f {x 2 +y 2 ) 
2 \ tanh(a)i) sinh(o;i) 



2xy 



(18) 



(19) 



is the well-known heat kernel, or the integral kernel of the operator exp(— tH), where 

Heat IVP-B. Consider the heat equation with quadratic potential and time variable coeffi- 
cients 



™ = l -e-^-^xH 
dt 2 dx 2 2 ' 

®{x,t)\ t =Q = <&(x, 0), -co < x < oo, 
where £i(t) = e 7 *, 7 > 0, and uj 2 {t) = — uj 2 , ui > 0. The corresponding IVP © is 

r + 7 r - w 2 r = 0, r(0) = r 7^ 0, r(0) = 0, 

which has solution 

r(t) = r o^e~^ cosh[Sli + /3], 



(20) 



(21) 



with = ■ x /^ 2 + (7 2 /4), and /3 = tanh _1 (7/2fi). Therefore, 



rj(x,t) 



c/i[m + /3], r(t) = ( tanh[Qt + /3] - 



2fi 



and according to (jlll) . the IVP ((20|) has general solution 



x exp 



2-kuj cosh[fti + P]r(t) 

1 / VLe^^x 



-e 7 * ( tanh[fk + /3] - X | .,- 



x / exp 

' — GO 



$(y,0)dy. 



2r(t) yuicoahpt + 0] 

In particular, when the initial function is $(2;, 0) = S(x), we have the fundamental solution 

Q 



(22) 



*(z,t) 



27ro) cosh[fti + 0\r(t) 



x exp 



fi 2 



=7* 



x exp 



1 



w 2 r(t) cosh 2 [fit + /3] 



or 



fie^/2 



2ttw cosh[fit + $]r(t) 



x exp 



l -e^ { u 2 r(t) + -„ ^ 

2 V uj 2 T (t) cosh 2 [tit + (3} 



.(23) 



Note that, when 7 — > 0, one has £1 — >• Q and /3->0, so that solution (1231) tends to solution (IT71) of 
the constant coefficient Heat IVP- A. The effect of the variable coefficients in IVP ((20)) . in particular 
the exponentially decreasing diffusion coefficient, can be seen by comparing FigJTJ-a, which shows 
solution (fT7|) and FigfTJ-b showing solution 
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Figure 1: a) Fundamental solution of Heat IVP-A at times t = 0.1, t = 1, t = 2, when u> = 1. 
b) Fundamental solution of Heat IVP-B at times t — 0.1, t = 1, £ = 2, when w = 1, 7 = 2. 



Another approach to solve the IVP © is using the Evolution Operator Method, also known 
as the Wei-Norman algebraic method, 27 , according to which, the evolution operator can be 
represented as an ordered product of exponential operators containing single generators of a Lie 
group. We obtain the following result. 
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Proposition 2.2 Solution of the IVP (0) is formally given by 

=W(t,toMx,t ), (24) 

where the evolution operator is 

p[t) 2 \ ( d , M _/V(t) d 2 



w(t,t ) = exp ( -^r x ) x ex P ( 9W x fa + 2H xcxp i 2 (25) 



and t/ie auxiliary functions are found in terms of solution r(t) of the IVP f3J) as 

p{t)=n{t)^§r; r(f)=r 2 (t ) f f r(t ) = 0; = - 1* 

Proof: The proof follows the same lines as for the Schrodinger equation in paper Ref . [24] . Indeed, 
for the IVP IJ5J, the evolution operator W(t,to) can be found by solving the operator equation 

^-W(t, t ) = f(t)W(t, t ), W(t , t ) = I, (27) 
at 

where 

f = m & + eSS T^ = «-W)*- - " (t) " J(t)Jf+) - <28) 

and 

K- = — — — , K+ = -x 2 , K = -(x— + -) 

2dx 2 ' + 2 ' ' 2 V dx 2' 

are operators which satisfy the commutation relations 

[k-,k+] = 2k Q , [k+,k ] = -k+, {k_,k } = k_, 

of the SU(1, 1) algebra. Assuming 

W(t,t ) = exp(ip(t)K + ) x exp(ir(£)X_) x exp(2g(t)k ), (29) 

and taking time-derivative of W, we get 

■^-WYi, t ) = \(ip - i2pg - ip 2 Te- 2g )K + + (iTe- 2g )k_ + (2g + 2pf e^ 9 )^] Wit, to). (30) 
at L J 

Comparing (I2T1) and (|30p . it can be seen that W(t, to) is a solution of (|27p. if the auxiliary functions 
p(t), r(t) and <?(£) satisfy the non-linear system of equations Q , © , ([5]) , which solution is given by 
(|26|). Finally, since IV satisfies ([27). it is well known that the function $(x,£) = W(t, t )$(x, t ) 
is a formal solution of the IVP ©. □ 

Combining the results of Proposition 12.11 and Proposition 12.21 it follows that the evolution 
operator W(t, to) given by (|25|) has integral kernel K t (x,y) formally given by ([T4|. i.e. 



W(t,t Q ) . = / K t (x,y) . aly. 

J — oo 

Moreover, using that 

^ *r(t) d 2 

exp[Ax— ]/(x) = /(e A x) and exp[-^— i ) = </?(a:, t(0)> 



where ip is solution of the classical heat IVP (|6]), we have 



$(z,t) 



W{t,t )${x,t ) 



exp 



P{t), 



d 1 



x 2 ) x exp ( a(t)(x— + 77) ) x exp 



-(t) d 2 
2 dx 2 



exp(ig(t)) x exp(-^ 2 ) x <p(e°Mx,T(t)), 



(31) 



and writing the functions p(t), r(t) and <?(t) as found in (|2l))) . one can see that solution (|3~Tj) 
obtained by the evolution operator method coincides with the previously found solution (U). This 
establishes the equivalence of both approaches. 

Proposition 12.21 can be easily applied when the initial function is given in the form of a power 
series, or in particular if it is a linear combinations of exponential functions such as cosh and sinh, 
using the well known relations, see [2"8] . 

d 2 , „ „ d 2 



exp[A— jx" 



H m (x, A), exp[X^]e x = eV, 



(32) 



where H m (x, A) are Kampe de Feriet polynomials defined by 



[m/2] 

H m (x, A) = ml 



fc=0 



fe!(m-2fc)!* 



iT m (z,0) = x" 



Therefore, using Proposition 12.21 exact solutions of the following parabolic IVP's are found. 
We remark that, these IVP's may not have direct physical meaning, but will be useful in Sec. 4, to 
obtain special exact solutions for the corresponding variable parametric Burgers type models. 

IVP 2.1. The IVP © with initial condition $(x,t ) = e~ cx cosh[Ax + c ], -00 < x < 00, 
and A, c, cq real constants, has solution 

' e (A — c)x+c 1 e -(A+c)x-c ' 



$(x,t) = W(t,to) 



<r(to) 



V »•(*) 

x cosh 



x exp 



2r(t) 1 



x exp 



,A 2 



)r(t) 



x exp 



r(tp) ' 
° r(t) . 



A(^x-cr(i))+co 



(33) 



where r(t) is solution of IVP (01, r(i) is as given in (|2"B|). 

IVP 2.2. The IVP © with initial condition $(x,t ) = x r ' 
has solution of the form 



$ m (x,t) = W(t,t )ai m = exp 



-00 < x < 00, m = 0,1,2,..., 



x ff m ( e 9(t) a;, £r(t) 



r(t) 



x exp 



2r(t) ' 



x ff„ 
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General IVP 2.3 The IVP © with general initial condition $(x,t ) = J2m=o a mX m , has 
formal solution 



${x,t) = ^a m $> m {x,t) 



m=0 



r(t) 



x exp 



2r(t) • 



r(t) ' 2 l j 



(35) 



where $ m (x, t) is given by (1341 . 



3 Variable parametric potential Burgers equation 

Now we consider the initial value problem for a variable parametric potential Burgers equation of 
the form 

dt n{t) 2 [ dx' 2fi(t) dx 2 2 ' (3 6 ) 
T(x,t)|i=i = T(a:,to) . 

We note that, if we write &(x,t) = exp(— fj,(t)Y(x, t)), then the variable parametric parabolic 
equation © transforms formally to the variable parametric potential Burgers equation (I36|) . On 
the other hand, writing ip(rj,T) — exp(— h(rj, r)) the standard heat equation transforms to the 
standard potential Burgers equation. Using this relations and Proposition 12.11 we formulate the 
following propositions. 

Proposition 3.1 The IVP h3b}) for the variable parametric potential Burgers equation has solution 
of the following forms: 

a) 

T <- " —m* ( w) W) hMx - '>• T{t)) <37) 

where r(t) is solution of IVP (0), 

v ( x ,t) = ^-x; r(t) = r 2 (t ) f f r(t o ) = 0, (38) 

r(t) J (C) 

and h(r),T) satisfies the IVP for the standard potential Burgers equation 

1 2 1 

kr + 2 " = 2 kr,v ' (39) 
h( V ,0) = fi(t )T( v ,t ) 



b) 



where tj,t are as defined in part (a), and ip(rj,r) satisfies the IVP for the heat equation 



(p(r), 0) = exp[-fi(t )T(j], t )}. □ 



(41) 
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Using the general solution (fit)]) of the heat equation and the relation h(rj, r) 
can easily write the analytic solution of IVP (131)1) , that is 



1 



exp 



-/j,(t )r(y,t ) 



2r 



dy 



In ip(rj, r), one 



(42) 



As a result, the general solution of IVP ((36)) is found in the form 
T(x,t) = 



1 ln fr(t 



2/x(t) 
1 



In 



r(t) 
1 



2r(t) a! 



exp 



-fj,(t )T(y,t ) 



(n(x,t) - yf 
2r(t) 



dy 



(43) 



y27rr(t) 

Next proposition gives direct relation of the variable parametric potential Burgers IVP with 
the variable parametric parabolic IVP. 

Proposition 3.2 The IVP LS6}) for the variable parametric potential Burgers equation has formal 
solution 

T(x,t) = 1 \n^(x,t), (44) 

fi(t) 

where <fr(x,t) satisfies the IVP for the variable parametric parabolic equation 

<9$ 1 <9 2 $ n(tW(t) 2 

dt 2fi(t) dx 2 2 ' ' (45) 

*(x,t ) = e3tp(-/i(to)T(a:,to)). □ 

4 Variable parametric Burgers equation 

Finally, we discuss the IVP for a one-dimensional variable parametric Burgers equation 



BU 



1 d 2 U 



dt ' n(t) U + U dx 2/x(t) dx 2 



\t)x , 



(46) 



U(x,t)\t=t = U(x,to) , — oo < x < oo 

where T(t) = fi(t)//j,(t) is the damping term, D(t) = l/2/i(i) is the diffusion coefficient, and 
F(x,t) = -u 2 (t)x is the forcing term which is linear in the space variable x. 

In [22], it was shown that solutions of the IVP (|4^| can be obtained in terms of solutions to 
the standard Burgers equation. For completeness and comparison, in Proposition 14.11 we outline 
this approach. Using it, special exact solutions such as generalized shock and multi-shock solitary 
waves, triangular waves, TV-waves and rational type solutions were found and discussed for forced 
Burgers equations with constant damping and exponentially decaying diffusion coefficients. 

Proposition 4.1 The IVP for the variable parametric Burgers equation J^6] ) has solution in the 
following forms: 



a) 

where r(t) is a solution of IVP 



u(x, t) = ^x + r }l a ) - v (n(x, t), r(t)) 



r(t) n(t)r(t) 



V (x,t) = r -^lx; r(t)=r 2 (t ) " 
r(t) 



r(to) = 0, 



(47) 



(48) 
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and the function V(r], r) satisfies the IVP for the standard Burgers equation 

Vr + VV v = \V m , 

V( V ,0)=^t o )U(n,t o ) 



(49) 



b) 



r(t) /i(t)r(t) (^(?7(x,t),T(t)) ' 
where 77, r are as defined in part (a), and ip(rj,T) satisfies the IVP for the heat equation 



(50) 



(51) 



'fir 2 ' 

pfa, 0) = exp [- p nfoMZ, ioK] . □ 

The well known solution (fTUf of the IVP (|5ip and the Cole-Hopf transformation V" = —(frj/ip, 
lead to solution of the IVP (|49l) for the Burgers equation 



V(r,,r) 



Coo (^) exp 




dy 


/•CO 

Loo eX P 


- 




dy 



Therefore, using Proposition 14.11 one can find formal solution of the IVP (|46)l for the variable 
parametric Burgers equation in terms of solution r(t) of the IVP ([3]), that is 



+ 



r(to) ' 


coo ( 1 /r(t ) 
J -00 \r(t) \ r(t) X 




exp 


/ 1 ^(to) 
V2r(t)U(t) 


s-y) a +/V(W(l/',toW)" 


dy 


tx(t)r{t)_ 


poo 

J -00 ex p 




(ar(t) 




h J v pL{t )U{y',t )d') 


dy 



Next proposition establishes direct relation between the variable parametric BE and the vari- 
able parametric parabolic equation via a generalized Cole-Hopf transformation. This gives us an 
alternative way of solving the IVP (|4l)| . 

Proposition 4.2 The IVP J46| ) for the variable parametric Burgers equation has solution of the 
forms: 

a) 



d 



U(x,t) = —T{x,t), 
where T(x,t) is solution of the IVP for the variable parametric potential BE 



1 d 2 r 



7 -(t) 



h i-^T + -( — r = 

dt n(t) 2 y dx' 2n(t) dx 2 
T(x,t ) = f x U(Z 1 t )dt. 



b) U(x t t) = --^£.Qn*(x,t)), 
where <f>(x,t) satisfies the IVP for the variable parametric parabolic equation 
( d<$> 1 <9 2 $ , /*(t)w a (t) 2$ 



dt 2/i(t) dx 2 2 
^ $(x,t )=exp(-/i(to)r^,*oK). □ 



(52) 



(53) 



(54) 



(55) 
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Proof can be done by direct calculation. In what follows, results of Sec. 2 are combined with 
Proposition 14.21 - (b) to obtain analytic and exact solutions of some Burgers IVPs. 
Burgers IVP-A. Given the Burgers problem with linear external potential 



1 9 

U t + UU x = -U xx + lj 2 x, 
U(x,t)\ t = = U(x,0), 



(56) 



-00 < X < 00, 



where Q > 0, it is not difficult to see that its linearization takes the form of the Heat IVP-A. 
Therefore, its solution is U(x,t) — —d x (lsi'^(x,t)), where ^(x,t) is given by (IT6"|) . Explicitly one 
has 



U(x,t) — Q tanh(u)t)i£ + sech(Cot)- 



roo ( sech(M)x-y \ \ ( sech(0it)x~y) 2 ry r,^ 

J-oo\ tanh^i)/^ J eX P [ ^ 2tanh(£i)/" + J U 



dy 



dy 



which is a known solution, Refs.|29].|3U]. 

Burgers IVP-B. The forced Burgers problem with damping term and exponentially decaying 
diffusion coefficient 

1 

s '"lA™. + u'x, Co > 0. 

(57) 



U t +-/U + UU X = -e^U xx + Co 2 x, Co > 0, 



U(x,t)\ t=0 = U(x,0), 



-oo < x < oo, 



is linearazible in the form of Heat IVP-B. Therefore, its solution is U(x,t) — — e '' t d x (ln'^(x,t)), 
where \&(x, t) is given by (|22D . Explicitly one has 



f7(x,i) = (--^ +fJtanh[f2i + /3]) a; 



a) cosh[fit + 0\ r(t) 



n e ^ 2 x 



Co cosh[fM + j3] 



V exp 



2r(t) \ Co coshfit + f3] 



n e ^/ 2 x 



$(y,0)dy 



exp 

I VP 4.1. The I VP (gBJ with initial condition 

1 



2r(t) \ wcosh[m + /3] 



$(y,0)dy 



l/(a;, to) = 



A*(*o) 

has single shock solitary type solution 

1 d 



(c — Atanh [Ax + cq]) , — oo < a; < oo, 



U(x,t) 



/j,(t) dx 



(ln$(x,t)) 



: — A tanh 



.4 



r(t) 



x ~ cr(t) + c 



r(t) /i(*)r(i) 

Here we used that &(x,t) is a solution of the parabolic IVP 2.1 given in Sec. 2. 
IVP 4.2. The IVP (gBJ) with initial conditions 

m 1 



(58) 
(59) 



U m (x,t ) 



M(*o) £ 



772 — 0, 1,2, ... 
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has rational type solutions 

1 d 



U m (x,t) 



r(t) 



m r(to) 



r(t) r(t) rr f rW 



m = 0,1,2, 



(60) 



where we used that $ m (x, t) is a solution of the parabolic IVP 2.2, and J^H m (x, A) = mH m —i(x, A). 
General IVP 4.3. The IVP (@6} with the general initial condition 

15 



E7(a;,t) 



M(*o) ^ 



n=0 



has formal solution given by 



^ 19 „ _ xx f(t) 1 d 

U(x,t) = __(]n$(a;,t)) = 



In a m H rr 



r -Mx l -r(t) 
r(t) ' 2 W 



/i(t) cte ' r(i) /i(t) 9x 

where we used that t) is a solution of the parabolic IVP 2.3. 

We note that the above special IVP's were discussed in [35], but using Proposition 14.11 Here, 
we illustrate the application of Proposition 14.21 (b). and show that both approaches lead to the 
same results. 



5 Summary 

In this work, variable parametric Burgers type models are introduced and two different approaches 
for solving the corresponding initial value problems are given. The first approach is transforming 
the variable parametric model to a standard constant coefficient model. The second approach is 
a direct linearization of the variable parametric Burgers model to a variable parametric parabolic 
model. Both approaches and the relations between them are summarized in the following scheme. 



Variable Parametric Models 



U t + ^§-U + UU X = ^7^U XX - oj 2 {t)x 
H{t) 2/4*) 



T, 



u(x,t) = r x (x,t) 
X 



\ x 

2n{t) xx 



T(x,t) 



I 

-^\n$(x,t) 
I 



** = 



2/x(t) 



Standard Models 







vv n 


= \v m 




I 




v(v, 


r) = 
I 


h v (r),T) 


h T + 






I 
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At a final stage, the problem of solving the variable parametric models introduced in this article, 
reduces to that of finding solution of a corresponding second order linear ODE with time dependent 
coefficients. This will allow us to study a wide class of exactly solvable Burgers type models related 
with the classical Sturm-Liouville problems for the orthogonal polynomials, or special functions. 
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